Surface tension and fingering of miscible interfaces 
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Abstract 

Experiments on miscible, buoyantly unstable reaction-diffusion fronts and non-reacting 
displacement fronts in Hele-Shaw cells show a fingering-type instability whose wavelengths (X ) 
are consistent with an interfacial tension (X) at the front caused by the change in chemical 
composition, even though the solutions are miscible in all proportions. In conjunction with the 
Saffman-Taylor model, the relation X ~ K lx, where x is the interface thickness and K = 4 ± 2 x 
10‘ 6 dyne, enables prediction of our measured values of X as well as results from prior 
experiments on miscible interfaces. These results indicate that even for miscible fluids, surface 
tension is generally a more significant factor than diffusion in interfacial stability and flow 
characteristics. 
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Introduction 


A Hele-Shaw cell is perhaps the simplest system in which multi-dimensional convection 
is present and is used extensively as a model hydrodynamic system [1], Buoyant convection and 
viscosity contrast cause instability of displacement fronts between non-reacting immiscible [1] 
and miscible [1,2] fluids as well as propagating reaction-diffusion fronts [3, 4, 5] in Hele-Shaw 
cells. A key instability characteristic is the wavelength (A,*) of the first fingers observed on an 
initially flat interface, given by 27t/k , where k is the wavenumber corresponding to the 
maximum growth rate (a*) of infinitesimal interfacial disturbances. To obtain a finite value of 
k\ a mechanism is needed to suppress disturbances of large k. This mechanism is thought to be 
surface tension for immiscible interfaces [1] and diffusion for miscible nonreacting [1, 2] or 
reacting [3, 4, 5] interfaces. 

Saffman and Taylor [6] derived a a-k relation for a non-reacting immiscible interface 
between a displaced fluid (subscript “1”) and a displacing fluid (subscript “2”) in a Hele-Shaw 
cell: 


<7 = 

MiZlh u J 

Pi ~ P2 

^gw 2 cos(0) 


Vi+th \ 


J 12 


S * 2 jJ 
12(a*i + fh ) 


( 1 ) 


where p is the dynamic viscosity, g the gravitational acceleration, p the density, 0 the angle of 
the cell from vertical, U the speed at which the displacing fluid is forced into the displaced fluid, 
£ the surface tension and w the cell gap. For the conditions studied in this work the temperature 
and concentration changes across the fronts are far too small to cause the viscosity contrast term 
to be competitive with the buoyancy term; in this case Eq. (1) predicts 
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where 5 = (p,— P 2 VP 1 is the dimensionless density contrast. 


For miscible fronts it is thought that surface tension does not plays a role and instead 
diffusion suppresses high-k disturbances [1]. For reaction-diffusion fronts the cr-k relation with 
Z= 0 and m - p 2 = 0 [4] is the same as Eq. 1 with U replaced by the front propagation speed (s) 
and the surface tension term replaced by k 2 D, where D is the reactant diffusivity. This leads to 
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where Ci = 9671 and C 2 = 1/2304. For reaction-diffusion fronts the interface thickness (x) D/s 
[3, 4] is time-independent. For non-reacting displacement fronts x increases with time, but a o-k 
relation has been derived assuming initially x = 0 [1,7]. This relation leads to Eq. (3) again but 
with Ci = 967t/(V5-2) and C 2 = 0.000391. 

In this Letter we show that even for miscible fronts, a weak effective surface tension 
exists that is generally a more significant mechanism in wavelength selection (via Eq. 2) than 
diffusional effects (Eq. 3) because it results in a much smaller a* than would occur were Z = 0. 
While prior experiments have noted evidence of surface tension at miscible interfaces [ 8 , 9, 10], 
a difficulty in quantifying this evidence exists because for displacement fronts x increases over 
time due to diffusion. In contrast, for reaction-diffusion fronts x is constant, thus a constant Z 
may exist. We will show that Z and thus instability characteristics can be predicted for both 
types of interfaces using a single simple model. 


Experimental apparatus and procedures 

Experiments were conducted in Hele-Shaw cells made of transparent Lexan. The cell 
dimensions were 200 mm x 200 mm x w, with 0.18 mm < w < 3 mm. The estimated Reynolds 
numbers = (<J*/k*)w/v of the resulting flows in the initial stage of instability is less than 4 for all 
cases studied, thus laminar Pouseille flow prevails. 

For reaction-diffusion fronts, the hydrosulfite-iodate autocatalytic system [11] was used, 
which has the overall reaction S 2 C> 4 2 ’ + IO3' + H 2 0 — » 2 SO 4 2 + I + 2 H . A typical composition 
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was [S 2 0 4 2 ] = [IO 3] = 0.015 M with NaOH added to obtain pH = 12, for which s = 0.020 cm/s. 
Values of s between 0.0015 and 0.042 cm/s were obtained by varying the chemical composition. 
The fronts were visualized using the natural color change or enhanced using pH indicators. Due 
to the slight temperature rise and composition changes across the front (discussed later), the 
products are less dense than the reactants, thus upward-propagating fronts are buoyantly 
unstable. The experiments were conducted by initiating fronts at the top of the cell and, after a 
flat downward propagating front developed, sealing the cell and inverting it. 

For non-reacting displacement fronts, the experiments were conducted in a manner 
similar to Wooding [2], The lower part of the cell was filled with the more-dense fluid, a 
horizontal barrier strip was inserted, and the upper part was filled with the less-dense fluid. The 
barrier was pulled horizontally out of the cell to start the experiment. 

Values of X were determined from images of the interface in the early stage of fingering, 
either by counting the fingers over a known width of the cell as in [2, 5], or determined from the 
peaks of spatial Fourier transforms of the fronts. There was no statistically significant difference 
between the two methods. The typical standard deviation of X * in one experiment was 0.2A.*, 
which exceeds the uncertainty in the measurement of X * itself. Measured standard deviations 
will be shown below for displacement fronts. 

For reaction-diffusion fronts, the (temporary) temperature-induced fractional density 
change 67 (typically 6 x 10 4 ) was determined by calculating the thermodynamic temperature rise 
and applying the thermal expansion coefficient of water. This rise was checked experimentally 
using a thermocouple immersed in the solution; the difference was less than 10% in all cases 
tested. The (permanent) concentration-induced fractional density change 8c was measured using 
a home-made dilatometer. 8r/6 c varied from about 20 for [IO3 '] « [S 2 0 4 2 ] to about 4 for [IO3 ] 
> [S2O4 2 ]. Except as noted, all calculations for reaction-diffusion fronts employ 8 = 8r + S c . 
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Results - reaction diffusion fronts 

Shortly after inversion the fronts began to exhibit fingers, a typical image of which is 
shown in Fig. 1 (upper). Figure 2 (left) shows X for varying Peclet number Pe = Uw/D, where 
U = a*/k* is the phase velocity of the most amplified linear mode according to Eq. 2 (o /k from 
Eq. 3 yields the same U and thus Pe except for a factor of 3/4) and D = 1.4 x 10 cm /s is the 
mass diffusivity of I0 3 ' [12]. With this definition Pe = (l/36)(pg8cos(0)w 3 /pD), which can also 
be interpreted as a Rayleigh number. Note that X* does not change systematically over more 
than two decades of Pe. The only experimental parameter that affected X substantially was 9, 
for which X* ~ l/- N /cos(0) (Fig. 2, right). For horizontal propagation the interface was smooth, 
verifying that the instability is driven by buoyancy. 

These results contrast markedly to the scalings of Eq. 3, where for a given composition 
(fixed p, 8, £, v, D), X* ~ w 2 /cos(0). Moreover, Eq. 3 predicts X*/w = (87t/3)Pe, which is lower 
than observed values by factors of 30 to 1500. This indicates that some process other than 
diffusion suppresses high-k modes. In contrast, our results are entirely consistent with the 
scalings of Eq. (2), i.e. X* ~ w°/Vcos(0) , suggesting that high-k modes are suppressed by 
surface tension — even though the reactants and products are miscible in all proportions. For a 
typical X* = 10 mm, 8 = 6 x 10 -4 and 0 = 0, Eq. (2) predicts £ » 5 x 10 3 dyne/cm, which is 1.4 x 
10 4 times smaller than that between water and air at 22°C. Thus, a rather weak surface tension 
could cause the observed behavior. 

Cahn and Hilliard [13] and Davis [14] suggest that for interfaces between miscible 
regular solutions having different solute concentrations 

(4) ' 

where K is a constant, C is the local solute concentration and AC is the total concentration change 
across the interface. (Surface tension could also result from the thermal gradients caused by heat 
release, but the thermal diffusivity a » D for most liquids, consequently the thermal T is much 
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larger than the compositional x, and thus would result in a much lower induced E.) For reasons 
not yet clear to us, our results and those of other investigations for both reaction-diffusion fronts 
and displacement fronts are considerably better described by neglecting the AC effect and 
presuming E = K/t, consequently we will employ this relation instead. 

For the aforementioned I = 5x 10' 3 dyne/cm and a typical s = 0.02 cm/s, we infer K » 
ED/s = 3.5 x 10' 6 dyne. To determine a best-fit K for all our data, the effects of heat loss to the 
cell walls must be considered as follows. It is readily shown that the total thermal energy (thus 
thermally buoyant fluid) in the system is proportional to the ratio (H) of the heat generation rate 
(=s 2 /a) to heat loss rate (=a/w 2 ), thus H = (sw/a) 2 . We expect that for H « 1, the thermal 
contribution to 5 is lost via heat transfer to the walls, thus 8 = 5c, whereas for H » 1, both 
thermal and compositional contributions are important, thus 8 = 8c+ 8r. This requires 
A, (experiment)/ A (theory, 8c+5t) be larger than unity (corresponding to an overestimate of 8) for 
small H, but decrease to unity at large H, and that ^‘(experiment)/ A*(theory, 5c) be close to unity 
at small H, but decrease to less than unity for large H. A best fit of our data to these 
requirements suggests K = 2.0 x 10' 6 dyne (Fig. 3). 

Results - displacement fronts 

Shortly after removing the barrier strip, the fronts started to exhibit fingers, a typical 
image of which is shown in Fig. 1 (lower). Figure 4 shows that A*/w was nearly constant (= 10) 
for a wide range of 5, w and 0 and different solutes, unlike reaction-diffusion fronts where A* 
itself was nearly constant for 0 = 0. The prediction based on the diffusion model (Eq. 3) is A*/w 
= 8n/(3(V5-2)Pe), which is far smaller than experimental values except at very small Pe, thus we 
again consult the surface tension model. For non-reacting interfaces, the model must be 
modified since x and thus E is not prescribed a priori by the balance between chemical reaction 
and diffusion. Since diffusion thickens the front while buoyant convection thins it, we estimate x 
as D/U where again U = o*/k\ and thus E* = Ko*/Dk*. Equation (1) then leads to 
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For our conditions (p. = 0.010 g/cm-s and D = 6.0 x 10 ^ cm /s for ICMn 04 [2]), using K — 2.0 x 
10 6 dyne as for reaction-diffusion fronts, Eq. 5 predicts X /w — 10.4. For Pe > 50, this is close to 
our measured values (Fig. 4) for a two-decade range of Pe, as well as a data point from Wooding 
[2], who mainly examined the nonlinear growth of large amplitude fingers but for two cases (0 = 
34. 1", 87.1°) reported an initial X*. 

Our apparatus did not permit low enough Pe to examine diffusion-dominated fingering, 
which should result in X*/w given by Eq. 3, but a transition to diffusion-q#ecte</ behavior is seen 
in Fig. 4 at Pe < 50. Our data are consistent with Wooding's data point in this regime. The 
transition is expected when values of a from Eqs. 3 and 5 are equal, which corresponds to Pe = 
^jpD/mClK ~ 43 for our conditions, which is close to the observed transition. The Pe « 1 
regime has recently been studied extensively [15] and good agreement with Eq. 3 is found. 

Comparison to prior works 

Carey et al. [5] studied rising reacting-diffusion fronts with very small s (0.00035 cm/s) 
and 8 (4.0 x 10‘ 5 ) in Hele-Shaw cells with 0.40 < w < 0.94 mm. They claim that their measured 
X* are consistent the diffusion model (Eq. 3), however, their Fig. 2 shows that for w > 0.5 mm, 
the data are fit equally well by a constant X* = 5 mm, which is consistent with the surface tension 
model. Using Eq. 2 for their test conditions then implies Z = 8.3 x 10' 5 dyne/cm and K = 3.3 x 
10‘ 6 dyne. The transition from diffusion to surface tension controlled behavior should occur 
when a* from Eqs. 2 and 3 are equal. This corresponds to Pe = (1024/27)(pD 2 /wKs) or w = 0.54 
mm for their conditions, which is close to their apparent transition from X* ~ w ' 2 to X ~ w°. 

Mayer and collaborators [16] measured an apparent £ at dissolving isobutyric acid/water 
(IBW) and cyclohexane/methanol (CM) interfaces using a capillary-wave technique. X varied as 
M^fDt ~ 1/t, where t is the time lapse since the mixture was heated from a two-phase 
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(immiscible) to a one-phase (miscible) equilibrium condition, which is consistent with the 
relation £ = K/x. For IBW at t = 5700 s, with a measured £ = 1.6 x 10' 4 dyne/cm and a reported 
D = 1.8 x 10 8 cm 2 /s, we infer K = 1.6 x 10 6 dyne. For CW at t = 2 x 10 4 s, £ = 8 x 10" 5 dyne/cm 
and D = 1.4 x 10' 8 cm 2 /s, we infer K = 4.2 x 10' 6 dyne. 

Matthiessen et al. [17] showed that for thin (= 1mm) horizontal layers of Belousov- 
Zhabotinskii reactants, the surface tension gradient at the free surface due to compositional 
changes across the front dominated buoyancy as a flow driving mechanism. They did not 
explicitly determine £, but found that the measured and calculated convection velocity profiles 
were best fit by a solutal Marangoni number of 0.3. While the £ gradient at the free surface is 
not directly related to £ across the chemical interface in the bulk fluid discussed in this work, 
their Marangoni number can be used to infer d£/dx = 0.0194 dyne/cm 2 . Using their presumed x 
= 0.018 cm, an apparent K * 6.3 x 10" 6 dyne can be inferred. 

Patterson [18] examined fingering in miscible fronts of water displacing glycerin in 
horizontal Hele-Shaw cells, where the viscosity contrast rather than buoyancy drives the 
instability. With pi » p 2 , g = 0 and £ = KU/D, where in this case U is a forced displacement 
speed, Eq. 1 predicts X / w = n-jKTju l D . With Patterson’s measured X*/w = 4, pi = 10.2 g/cm-s 
and D = 1 .7 x 10‘ 7 cm 2 /s [8], we infer K = 2.8 x 10 6 dyne. 

Discussion 

Our results show that interfaces between miscible fluids possess a weak surface tension 
that affects the interfacial stability and resulting flow, typically more significantly than diffusion. 
An empirical constant K = £x = 4 ± 2 x 10‘ 6 dyne describes a broad range of experimental 
results. This magnitude of K might be anticipated since the maximum attractive force between 
water molecules based on the Stockmayer potential is (504/1 69)(7/26) l/6 (5.2 x 10‘ 14 erg)/(2.6 x 
10’ 8 cm) = 4.8 x 10’ 6 dyne. Davis [14] recommended K = 2 x 10‘ 6 dyne for hydrocarbons but 
gave no justification. While we have not determined the limits of applicability of our K value, 
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even for an immiscible water-air interface with x ~ 10' 7 cm, the prediction L = 40 dyne/cm has 
the right order of magnitude. 

We thank Drs. H. T. Davis, S. Davis, J. Goddard, T. Maxworthy and P. Petitjeans for 
helpful discussions. This work was supported by NASA grants NAG3-1523 and NAG3-2124. 
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Figure 1. Contrast-enhanced images of interfaces in Hele-Shaw cells. Upper: upward 
propagating reaction-diffusion front with s = 0.0038 cm/s, w = 1.0 mm, 8 = 0.00032, Pe = 63, 
field of view 5.0 cm wide. Lower: non-reacting displacement front of water over water/ethanol 
solution with KMnC >4 dye, w = 0.80 mm, 8 = -0.00036, 9 = 180°, Pe = 85, field of view 5.4 cm 
wide. 
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Figure 2. Observed instability wavelengths (X,*) for upward propagating reaction-diffusion fronts 
in Hele-Shaw cells. Left: as a function of the Peclet number = (l/36)(pg8cos(0)w 3 /pD), for 
vertical cell orientation (0 = 0). Mean and standard deviation for all points shown are 9.5 mm 
and 2.2 mm, respectively. Right: as a function of 0 for s = 0.02 cm/s and w = 1.0 mm. 
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H = (sw/a ) 2 


Figure 3. Comparison of measured and predicted (Eq. 2, using D = Ks/D, with K = 2 x 10' 6 
dyne) instability wavelengths for dimensionless density contrasts due to both thermal and 
composition changes (8t + 8c) and composition changes (8c) only, as a function of the heat loss 
parameter H. 
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Figure 4. Measured initial instability wavelengths (X*), referenced to the cell thickness (w), for 
non-reacting solute-water solutions over or under water, as a function of the Peclet number, 
along with the predictions of Eqs. 3 and 5 for K = 2 x 10 dyne. For the data set of varying 8, 
the solute is KMn0 4 , w = 1.3 mm and 9 = 0. For the data set of varying 9, 8, the solute is 
KMn0 4 , w = 0.80 mm and 8 = 0.0017. For the NaCl and NaI0 3 data points, w = 0.80 mm, 8 = 
0.0017 and 9 = 0. For the ethanol data point, w = 0.80 mm, 8 = -0.00036 and 9 = 180°. Where 
shown, error bars correspond to one standard deviation of the measured X for one test. 
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